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We show using bosonization techniques, that strong forward scattering interactions between 
one dimensional spinless Luttinger liquids (LL) can stabilize a phase where charge-density wave, 
superconducting and transverse single particle hopping perturbations are irrelevant. This new phase 
retains its LL like properties in the directions of the chains, but with relations between exponents 
modified by the transverse interactions, whereas, it is a perfect insulator in the transverse direction. 
The mechanism that stabilizes this phase are strong transverse charge density wave fluctuations at 
incommensurate wavevector, which frustrates crystal formation by preventing lock-in of the in-chain 
density waves. 



Interacting fermions in one dimension can exhibit Lut- 
tinger liquid behaviour where, in contrast to Fermi liq- 
uids, all the low lying excitations are collective modes 
0. The search for such non- Fermi liquid behaviour in 
higher dimensions prompted several authors to study the 
problem of coupled Luttinger liquids. But renormal- 
ization group (RG) studies, such as || found the cou- 
pling to destabilize the Luttinger liquid behaviour. It 
was however argued by Anderson et al. that, in spite of 
its relevance in the RG sense, electron transverse hop- 
ping may be incoherent allowing for a non- Fermi 
liquid in dimensions greater than one. In addition to 
these theoretical motivations, interacting Luttinger liq- 
uids could well arise in experimental systems like quasi 
one dimensional organic conductors and in "ropes" of 
nanotubes. Recently, coupled one dimensional systems 
have also emerged in the stripe phases of Quantum Hall 
systems in higher Landau levels || , and in the cuprates 
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In this paper, we will revisit the problem of coupled 
Luttinger liquids, and address the question of the exis- 
tence of a stable two dimensional phase that retains some 
of the properties of the one dimensional Luttinger liquid. 
We consider this issue in detail within the RG framework 
using bosonization, for the case of spinless fermions. Al- 
though most of the physical systems of interest are of 
spinful (rather than spin polarized or spinless) fermions, 
that case is technically more involved due to the presence 
of exchange interactions between the chains, and we post- 
pone its study to || . The case of an infinite set of coupled 
spinless Luttinger liquids is simpler, but should still cap- 
ture some of the general physics present in the spinful 
case. There is also a possibility that such systems may 
be directly realized by spin polarizing a layered quasi- 
one-dimensional system of low electron density with an 
in plane magnetic field; or in the spin polarized Luttinger 
liquid formed by Zeeman split quasiparticles bound to the 
superconductor vortex core j?],|8| . The new ingredient in 
our study is the inclusion of strong forward scattering 
interactions between the chains. We then find that an 
anisotropic phase, which is Luttinger liquid like along the 
chains but a perfect charge insulator in the transverse di- 



rection, is stable against a range of instabilities including 
those that could lead to a superconductor, crystal or 2D 
fermi liquid phase. Hence this new phase is a non-Fermi 
liquid phase in two dimensions that is highly anisotropic. 
In the direction of the chains the correlation functions 
have power law forms with nontrivial exponents, as in a 
Luttinger liquid. However, the relations between expo- 
nents are modified, as compared with completely decou- 
pled Luttinger liquids, by the strong forward scattering 
interactions between the chains. We will call this phase 
the sliding Luttinger liquid (SLL). It is analogous to the 
sliding phase found by O'Hern et al. |J in the related 
problem of XY models coupled by suitable gradient in- 
teractions, which motivated our approach. In contrast 
to p0[ , we consider transverse hopping operators when 
establishing the stability of the SLL. 

The stable sliding Luttinger liquid fixed points that 
we find occur close to an instability towards transverse 
charge density wave ordering. The instability occurs 
when the stiffness of the density fluctuation mode at a 
particular transverse wave- vector vanishes. We propose a 
physical mechanism that is responsible for stabilizing this 
phase - i.e. that the proximity to the transverse charge 
density wave state induces strong fluctuations of the lo- 
cal in-chain density, which frustrates crystal formation. 
This mechanism is most effective if the wavelength of the 
transverse charge density wave phase is incommensurate 
with the spacing between the chains. 

The existence of this phase may at first sight contradict 
previous results on coupled Luttinger liquids. However, 
previous approaches focused either on the simpler case 
of two interacting chains or on the perturbative ef- 
fect of these transverse interactions. Indeed, as we will 
show this phase can exist when strong forward scattering 
couplings between nearest and (at least) second nearest 
neighbour chains are taken into account. 

The sliding Luttinger liquid. We consider an 
anisotropic 2D system composed of parallel chains (la- 
belled by an integer i) with spinless Luttinger liquids in 
each chain. The bosonized form of the fermion operators 
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near the right and left fermi points are: 
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^R/L,»(a;) 



A/L,i(^) 



where e is some intra-chain cut-off, the Xr/M are the 
Klein factors and x is the coordinate along the chain. 
The Luttinger liquid on each chain is described by the 
usual lagrangian: 



(i) 
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[(d t Qi(x,t)) 2 - (d x Qi(x,t)f 
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where 0j = (<^£,i — <j>R,i)/\4ft, and -Ko is an interaction 
dependent constant with iio > 1 for repulsive intra-chain 
interactions. The sound velocity on each chain has been 
set to unity. In the following, we will also use the dual 
field $i(£c, t) = {4> L , % + (f>R,i)/V4n. 

We now add forward scattering interactions between 
the chains which correspond to couplings between the 
long wavelength components of the densities p%(x, t), and 
of the currents Ji(x, t): 

C mt = 2vr ]T [JiKfi J j ~ PiK^Pj] (2) 

where the interactions are assumed local in x and time. 
Using the bosonization relations \f^npi{x, t) — d x @i(x, t) 
and y/4nJi(x,t) — —dt&i(x,t), we obtain the bosonized 

form of the lagrangian C to t = J2i + ^ int f° r tne m " 
teracting Luttinger liquids : 

id 
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where the coupling matrices are defined by K(j — KoSij 4 
K/j,K?j = KoSij + K?j. By introducing Fourier trans- 
forms in the direction transverse to the chains, this la- 
grangian can be rewritten as 



C-tnt — 
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v(q±) 



\d t e q 



v(q±)\d x Q q 
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with the notation f = f n -S^- (the transverse chain 
spacing has been set to unity). The stiffness K(q±) is 
defined by K(q±) = y/K J (q±_)KP(q±) and the veloc- 
ity v(qx) = y/KP(q±)/K J (q±). Note that Lorentz in- 
variance corresponds to K p (q±) cx K J (q±) (like in the 
isotropic model studied in H). In that case, all modes 
have the same velocity. 

The lagrangian ([|) is invariant under the transforma- 
tions $j — > $i + Cj and Qi — > 0i + di where Ci and di 
are constants on each chain. We will thus call the cor- 
responding fixed point a sliding Luttinger liquid (SLL) 



fixed point pi| . From this symmetry we deduce that 
the total numbers of left (right) moving fermions on each 
chain are good quantum numbers and expectation val- 
ues of operators that change these - such as {ip' L ^L,j) 
for i ^ j - are necessarily zero in this phase. This 
corresponds to a perfect charge insulator in the trans- 
verse direction. Density (and current) correlations in the 
transverse direction are however nontrivial. For short 
ranged density and current interactions, they decay ex- 
ponentially with separation between the chains. The low 
energy modes are density oscillations (sound) with dis- 
persion E{q«,q±) = v(q±)\qu\ (where q\\ is the wavevector 
along the chain) which in general will propagate both par- 
allel and perpendicular to the chains. These modes can, 
for instance, transport heat perpendicular to the chains 
although the system is a perfect charge insulator in that 
direction. Correlation functions along the chains exhibit 
power law behaviour as in the Luttinger liquid. However 
the exponents now depend on the function K(q±_) rather 
than on a single number as in the case of completely 
decoupled Luttinger liquids (|l|). Therefore, relations be- 
tween exponents that are valid for decoupled Luttinger 
liquids no longer hold in this case. 

As this phase is described by a gaussian lagrangian, we 
can study perturbatively the possible relevance of various 
operators to ascertain its stability. 

Transverse hopping operators. In the usual case of 
two chains, it is easily shown that the most relevant op- 
erators correspond either to single particle (SP), particle- 
hole (CDW) or pair hopping (SC) 0. It is thus natural 
to first focus on these operators in our stability analysis. 
They are defined respectively by 

<5£ SP = £ % (<^aj + + h.c) (5a) 

<5£cdw = Y fcow (^LiV'M^jV'Rj + h-cj (5b) 
SC SC = ]T gg (vL<^h,Aj + >'■<■■) ( 5c ) 

ij 

Upon bosonization, the corresponding operators read: 

0% = cos V5r($< - *,•) cos y/^(Qi - %) (6a) 
OS JW = cosV4i(e i -e J ) (6b) 
Og = cos V^tt ($i - $j) (6c) 

The dimensions of these operators at the SLL fixed 
point (W) are readily evaluated and are given by : 



(N) _ I 1 ~ cos(Nq ± ) 

ICDW — ^ 



K(q±) 



(7a) 



via } = 2 / (1 - cos(N qi _))K(qs_) (7b) 

Jq ± 



(7c) 
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where JV = \i — j\. A necessary condition for this sliding 
LL phase to be stable is thus that 

&,r$ ) >2 and r,™+r$ ) >8 (8) 
for all JV > 1. 

Stability of the SLL phase. Before discussing the sta- 
bility of SLL fixed points, let us first discuss the do- 
main of validity of our RG approach. We are considering 
the scaling dimension of perturbing operators around the 
fixed points defined by the action (||). This action is de- 
fined provided the stiffness K(q±) is positive everywhere 
in [—it, 7r] . Our study is thus restricted to the correspond- 
ing subspace of K(q±). At the boundary of this subspace 
K (q±) vanishes for some [0, tt] Jl^l- The density cor- 
relations of transverse wavevector q^_ then diverge, sig- 
nalling an instability towards transverse charge density 
wave ordering. In the transverse CDW, (p q <^) ^ 0, so 
then the total charge density on each chain is a func- 
tion of the transverse position i |l3| ] . That the boundary 
corresponding to this transition will play a crucial role 
in the following analysis can be seen by inspecting the 
dimensions (Q) : ?7cr>w can be significantly increased by 
the presence of the pole in the integrand (root of K(q±)) 
with rjic^ being not much affected Hence it may be 
possible to have these operators irrelevant (i.e. with di- 
mension greater then 2), for parameters in the vicinity of 
this boundary. 
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FIG. 1. SLL fixed points for the model in (|). The al- 
lowed parameter range for which K(q±) is positive is within 
the curve ACBD. For each point shown, there exists a finite 
range of Ko for which the SLL is found to be stable against (a) 
CDW, SC and SP perturbing operators of eq.(|) with N < 4 
shown by the small dots and (b) including general four point 
operators with JV < 10 shown as squares. 

Stability analysis of a model K(q±) Here we shall con- 
sider the stability analysis for a concrete model of K(q±). 
We thus look for a natural restriction to a finite number 
of terms of the Fourier expansion, that may allow a stable 



SLL fixed point. The simplest case turns out to be 

K(q±) = K [l + Ai cos(g ± ) + A 2 cos(2<? ± )] (9) 

(As we will discuss later, the model with just Ai in- 
cluded does not possess a stable SLL phase). The re- 
quirement of positivity of K(q±) restricts the range of Ai 
and A2, as shown in the Figure [l], to the region within 
the boundary ADBC (henceforth simply denoted by B). 
The stability of the SLL to the perturbations (||) with 
1 < (JV = \i - j\) < 4 is first determined Q. The expo- 
nents in eq. (Q) are numerically evaluated and the results 
are shown in Figure ^ which is a two dimensional repre- 
sentation of the (Kq,Xi, A2) space. The points marked 
are those for which there exists a range of Ko values 
where the above operators are all irrelevant. The corre- 
sponding values of Kq are all greater than one (repulsive 
interactions on the chains). All the stable fixed points 
are thus found close to the boundary B. 

In an attempt to define stable RG fixed points, one may 
worry about more general perturbing operators. Indeed, 
inclusion of the operators (||) with 5 < JV < 10 does not 
significantly change the results. However, it turns out 
that general four point operators such as, for example, 
simultaneous hopping of fermions from chain i to i + 1 
and from j to j + 1 [ p!5[ , has a dramatic effect on the 
stability region which is now much reduced. Thus, these 
operators are found to be relevant at a large number of 
the previously stable points. There are however some re- 
maining fixed points, clustered close to the boundary B 7 
which are shown as squares in Figure |l|. 




FIG. 2. Range of Ko corresponding to stable SLL fixed 
points as a function of the transverse CDW wavelength 2-K/q^. 



Inset shows actual values of Ko 

Thus, as expected from the argument given above, all 
stable fixed points are found close to the boundary B. 
This can be physically understood as follows. 

Let us recall that anywhere on B the system is on the 
verge of a transverse CDW instability. Such a transverse 
CDW would frustrate crystallization of the fermions since 
(2/«f) _1 , which controls the spacing of the particles on 
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each chain, is now a function of the chain index i. Strong 
fluctuations of this kind prevent the locking in of density 
modulations along the chain, hence defeating the crys- 
tal instability and stabilizing the SLL. This mechanism 
will be less effective if the wavevector g° is commensu- 
rate with the transverse lattice i.e. if q\ = 2nm/n for 
some integers m,n. Then, the longitudinal density waves 
on chains separated by n can lock in and give rise to the 
crystal instability (provided n is not too large). This idea 
is confirmed on inspecting the range of Kq for which the 
SLL phases exist. In Figure [2] the corresponding range 
of Kq is plotted as a function of the transverse CDW 
instability period 2Tr/q^_ (that it is nearest to). Drastic 
reduction of the SLL stability occurs at commensurate 
transverse wavevectors. This also allows us to under- 
stand the absence of stable SLL in the model with A2 = 0. 
Then, the only existing transverse CDW has wavevector 
7r (for Ai = 1) and hence for strong enough repulsive in- 
teractions, we expect the longitudinal density waves on 
next nearest neighbor chains (N=2) to lock and lead to 
a crystalline instability. This turns out to be a correct 
expectation as, by including the N=2 CDW operator, we 
do not find a stable SLL in this model. 

The physical mechanism identified here allows us to 
generalize our results to models beyond the simple ones 
considered so far, including Luttinger liquids coupled in 
three dimensions. In general we expect the SLL to be 
stabilized for renormalized couplings in the vicinity of an 
incommensurate transverse CDW transition. 

An approach very similar to the one taken in this paper 
can be applied to bosonic models at any filling, coupled 
anisotropic spin chains, or vortex lattices in anisotropic 
superconductors. These, together with the more compli- 
cated case of spinful fermions will be discussed in a forth- 
coming publication ||. As disorder in known to modify 
significantly the behaviour of a single Luttinger liquid, it 
may also be highly interesting to study its effect on the 
new phase described in this paper. 

In conclusion we have studied in detail the occurrence 
of the sliding Luttinger liquid phase in a simple three pa- 
rameter model. We find that the region of stability of this 
phase is restricted to be close to the boundary where the 
harmonic boson theory breaks down and a soft mode ap- 
pears, that signals a transverse CDW instability. Strong 
incommensurate transverse CDWs are particularly effec- 
tive at frustrating crystal formation and enhancing the 
stability of the SLL. We expect these conclusions to apply 
also to more general models than the simple one consid- 
ered in this paper. 

During the completion of this work, we became aware 
of a similar, but nevertheless different, study by Emery 
et al. ftfity. 
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